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ABSTRACT 

We explore the possibility of light and superlight sterile neutrinos in the recently proposed Minimal 
Radiative Inverse Seesaw extension of the Standard Model for neutrino masses, in which all existing 
neutrino data can be explained. In particular, we discuss two benchmark scenarios with one of the 
three sterile neutrino states in the keV-range, having very small mixing with the active neutrinos to 
account for the Dark Matter in the Universe, while (i) the other two light sterile neutrino states could 
be in the eV-range, possessing a nonzero mixing with the active states as required to explain the 
LSND+MiniBooNE+reactor neutrino data, or (ii) one of the light sterile states is in the eV-range, 
whereas the second one could be superlight and almost mass-degenerate with the solar neutrinos. 
Such superlight sterile neutrinos could give rise to potentially observable effects in future neutrino 
oscillation experiments and may also offer a possible explanation for the extra radiation observed 
in the Universe. 



I. INTRODUCTION 



A large number of solar, atmospheric, reactor and accelerator neutrino experiments provide solid pieces of evidence 
for the neutrino oscillation phenomenon, and hence for nonzero neutrino masses and mixing [T]. Most of these 
experimental results can be understood by oscillations of the three so-called active neutrinos, which appear in charged- 
and neutral-current weak interactions in the Standard Model (SM). A global fit 013] to the three-neutrino oscillation 
data gives their best-fit mass-squared differences of Am 2 ol = 7.5 x 10~ 5 eV 2 and Am 2 tm — 2.4 x 1CU 3 eV 2 . However, 
there are a few experimental results, most notably the LSND anomaly [3] and more recently the MiniBooNE results [5], 
as well as the reactor anti-neutrino anomaly [5] , which cannot be explained within this three- neutrino-mixing paradigm 
and seem to require additional neutrino(s) with mass-squared differences of order 1 eV 2 . Such additional neutrinos 
cannot couple directly to the SM Z-boson, due to the LEP constraint on its invisible decay width [7] which allows 
only three active light neutrinos. Hence, these neutrinos must be sterile neutrinos [5], i.e. SM gauge singlets, and can 
only participate in the SM weak interactions through their mixing with the active neutrinos. 

In the light of recent predictions of a slightly higher reactor anti- neutrino flux [5], a global analysis of the short 
baseline (SBL) neutrino oscillation data favours the existence of more than one light sterile neutrino in the eV mass 
range [lOHM] . On the other hand, cosmological data also indicate a weak preference for additional light degrees of 
freedom: N c g = 4.34±0.87 [15] . which could be interpreted as another evidence for sterile neutrinos in the eV/sub-eV 
mass range [TH, 116) . 

The existence of light sterile neutrinos requires further theoretical justification, since unlike the active neutrinos, 
the sterile neutrino masses are not protected by the SM gauge symmetry. There are a number of interesting proposals 
in the literature jSj 1 to explain the lightness of the sterile neutrinos in various extensions of the SM, for nonzero active 
neutrino masses vn, (with i = e, /j,,t). In the simplest extension which realizes the type-I seesaw mechanism [TH], ur 
right-handed neutrinos v a R (with a = 1, 2, tir) are added to the SM some of which could in principle play the role 
of light sterile neutrinos, provided their masses are in the eV range [19] . The masses and mixing of the active and 
sterile neutrinos can be deduced from the (3 + hr) x (3 + hr) complex symmetric neutrino mass matrix in the flavour 



1 Light sterile neutrinos can also occur in theories of quantum gravity due to gravitational interactions involving global anomalies related 
to quantum torsion 1171 . 
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basis {{v iL ) c \v aR }: 

- ( Ml Z ) ■ <» 

where Md is the Dirac neutrino mass matrix, and Mr is the (B — L)-breaking Majorana mass matrix of the singlet 
neutrinos. In the usual seesaw approximation ||£|| = HM^M^ 1 ! <C 1, with ||£|| = ^/Tr(£t£) being the norm of the 
matrix £, the light neutrino masses are determined from the mass matrix 

M VL ~ -M D M^Ml , (2) 

with mass eigenvalues whereas the heavier masses are roughly the eigenvalues of Mn: mi^...,^+ nR 3> 01^2,3. 

In this case, the active-sterile mixing, given by O ~ {mi, 2,3/^4,5,. ..^J 1 ' 2 , is very small, thus limiting our ability 
to observe the singlet neutrinos, unless Mn £ 10 eV which are mostly excluded over a wide range by the present 
data [20]. 

It is therefore interesting to examine a theoretical framework yielding simultaneously light sterile neutrinos as 
favoured by the SBL neutrino oscillation data and a sizable active-sterile mixing to be observable in the current or 
near future experiments. An attractive possibility is the so-called inverse seesaw model |21j . where in addition to the 
right-handed neutrinos as in the type-I seesaw model, another set of SM singlet fermions S p l (with p = 1, 2, ns) is 
introduced. In this case, the mass matrix in M gets extended to the following (3 + tir + ns) x (3 + n R + ns) general 
mass matrix in the basis {{vir,) c , v a R, (S p l)j- 

( M D \ 

M v = Ml u R Ml , (3) 
\ M N ps ) 

where in addition to the usual Dirac mass matrix Md and the Majorana mass matrix ur in the v^-vr sector, we 
have added another Dirac mass matrix Mn and Majorana mass matrix fig in the vr-S singlet sector. Observe 
that the standard inverse seesaw model discussed originally in |21j is recovered, once we set the right-handed neutrino 
Majorana mass u R = in In this case, for \\fis\\ "C ||Mjv||, ||Md||, it is possible to have very light sterile neutrinos, 
e.g. in the presence of a /1 — r symmetry [22] or in theories with warped extra dimensions [23] . 

Another interesting realization of inverse seesaw models arises, when \j,r ^ 0, but ps = 2 - In this case, the rank of 
the mass matrix given by ^ reduces to 3 + ur, and the light neutrinos are exactly massless at the tree level. However, 
they acquire a small mass at the one-loop level which is directly proportional to the Majorana mass matrix fiR |25j 
and results from well-known SM radiative corrections involving the Z- and Higgs bosons [26]. This scenario, called 
the Minimal Radiative Inverse Seesaw Model (MRISM) in [53] , is very economical, as it does not require the existence 
of other non-standard scalar or gauge fields or other fermionic matter beyond the singlet neutrinos {v a R, S p l}- 

In this paper, we study the possibility of having light sterile neutrinos in the MRISM with three pairs of singlet 
neutrinos {viL^SpL} (for i,p = 1,2,3). In the limit \\/j,r\\ 3> ||Md||, ||Mjv||, after integrating out the three heavy 
singlet states fi, 2,3/1 with masses of order \\(j.r\\, we are left with six light mass states. Of these six light states, three 
should describe the mostly active neutrinos and the remaining three will be mostly sterile states. In particular, we 
discuss two benchmark scenarios. In the first scenario, all the sterile neutrinos are heavier than the active ones, with 
two of them having mass in the eV-range (as in the 3+2 models) to explain the LSND+MiniBooNE+reactor data, 
while the remaining one is in the keV-range to account for the Dark Matter (DM) in the Universe 27J. In the second 
scenario, one of the sterile neutrinos is in the keV-range as in the first scenario, and the second one is in the eV-range 
(as in the 3+1 models), while the third sterile neutrino is superlight and almost mass-degenerate with the electron 
neutrino v e L- We present numerical estimates of the model parameters for both of these benchmark scenarios. Finally, 
we also comment on another possibility discussed in the literature when all the three sterile states are in the eV-range 
(the 3+3 models). 



2 This structure may arise, for instance, in models where the Majorana masses for the S-fields are forbidden due to specific flavour 
symmetries |24| . 
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The paper is organized as follows. In Section [TTJ, we briefly review the MRISM and present the effective neutrino 
mass matrix after integrating out the heavy sterile states. In Section |III| we analyze the simple case of one single 
flavour, for illustration purposes. In Section |IV| we generalize our analysis to the three flavour case and discuss the 
experimental constraints on the mixing parameters. In Section [Vj we present our results for the benchmark scenarios 
mentioned above. Our conclusions are given in Section |VI| In Appendix A, we give a general parametrization of the 
unitary matrix used in our analysis. 



II. STERILE NEUTRINOS IN THE MRISM 



In the limit \\hr\\ ^> ||-Md||, ||Mjv||, the ^-fields decouple below the mass scale (1r, resulting in an effective theory 
with six neutrino states: v e _ f _ l ^ T L and £1,2, 3L- At the tree level, the effective neutrino mass matrix in the weak basis 
{vr,h,tL, Si,2,sl} becomes 



tree _ ( M Dl i£Ml M D ^~ R l Ml 



N _ 



M N 



Mfl 1 ( Ml Mj 



N 



(4) 



Note that one of the block eigenvalues must vanish, since the effective mass matrix has rank 3. The matrix 



■MJg can be block-diagonalized by a unitary transformation: 



0, 



3 

.-1 n/rT 



3 M^M' 

where the unitary matrix V has an exact representation in terms of an arbitrary matrix £ 



V = 



(i 3 + ec T r 1/2 c*(i3 + c T c*r 1/2 
-c T (i 3 +cc T )- 1/2 (i 3 + c T cr 1/2 



From ([5]), it is straightforward to obtain 

C = M D M N 1 and M = (1 3 + C f C) 1/2 M N . 



(5) 



(6) 



(7) 



At the one-loop level, the mass matrix given by Q receives an electroweak radiative correction proportional to 
US], which in the limit ||/Xr|| 3> ||Md||, ||-Mjv|| is given by 



M 



1 — loop 

iff 



3 



aw 



167rm^ 



Mi? 



In 



/4 



t J 'R- m H 1 s \ m H/ VR-m^U 



3m| 



In 



Ml 3 



( M D y, R 1 XRf{x R )Ml 0s\ = ( AAf 3 

3 3 



O3 3 

where the one-loop function /(xr) is defined as 



f{x R ) 



aw_ 
16tt 



x H 



xr - x H 



In 



xr 
x H 



3x2 



xr - x z 



In 



xr\ 

Xz) 



(8) 



(9) 



with xr = pi R /m w , xh — ritfj/myy, %z = m'z/ rn wi assuming \ir = \xr\-z for simplicity. The functions /(xr) and 
x Rf( x R) are plotted in Fig. flj 

Thus, the full neutrino mass matrix in the basis {v e ,n,TL, S±_2,3l} is given by 



Meff = M^ c + M 



1— loop 
off 



M D ^ ( 1 3 + xrJ(xr)) Ml Md^M] 



-1 A.fT 
N 

M n1 jlr 1 MI M N fi R 1 Mj f 



(10) 
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FIG. 1: The analytic dependence of the one-loop function f(xa) and xnf{xa) on the Majorana mass scale fin. 



Since the loop- function Xnf(xjA is much smaller than unity for a wide range of (xr as shown in Fig. [T we can use 
the unitary matrix V given by ^ to perturbatively block-diagonalize the full mass matrix in ( |10[ ), as fo lows: 



V T M cS V 



(i 3 + CC t )" 1/2 AM(i + cC T )" 1/2 
[i 3 + C f C) 1/2 C t AAf(i + c*C T ) 1/2 

mil 

m T l2 M^M T + 77122 



(l3 + CC f ) 

J 1 M T + (la + Ck) 



" 1/2 AM(* (1 + C T C) _1/2 



-1/2 



C f AMC (1 + 



-1/2 



(11) 



Since ||mi 2 ||, HTO22II < 
eigenvalues: 



jM/j^M 1 



we can make the usual seesaw-like approximation to obtain the block mass 



Mi 
M 2 



~ 1 M r ) 



(12) 
(13) 



It is clear from ( 12 ) that the light neutrino mass block eigenvalue is proportional to the loop-correction factor AM, 
and vanishes as AM — > 0. Note that the second term in (12) represents an effective two-loop effect of 0((AM) 2 ), 
and next to the first term, it can be safely ignored. 



In order to obtain the active and sterile components for each mass eigenvalue, the mass matrix in (11) needs to 
be further diagonalized into a fully mass-diagonal form. Before doing so for the general three-flavour scenario in 
Section |IV[ we will first consider a simple one-flavour scenario for illustration purposes, which would help us to gain 
valuable insight. 



III. THE SINGLE FLAVOUR CASE 



In this case, we replace all the matrices in (12) and (13) by complex numbers. Then, the mass eigenvalues are 
simply given by 



777,1 



XRf{x R )m 2 N ( C 



i + ICI ; 



m 2 = ^(l + |C| 2 ) 



(14) 



while the mass eigenvectors are given by 



ni 

"2 



VS 



(i + KlW^f£ 



(15) 
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FIG. 2: The active and sterile components in the lighter neutrino state n% defined in ( 15 1 for the single flavour case 
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FIG. 3: The light neutrino mass eigenvalues mi and m,2 for the single flavour case with two different choices of xr = /ijj/mf^ 
and with |rrijv| = 10 _4 /xh (solid), 10 -5 fin (dashed). The horizontal red (dotted) line represents the tentative upper limit on 
the active neutrino mass scale. 



From (15), we note that in the limit ||£|| "C 1, the massless state is mostly the active neutrino vl, whereas in the limit 
HCII ^> 1, the massless state is mostly the sterile neutrino Sl- This generic property is illustrated in Fig. [2] 

The magnitudes of these eigenvalues are shown in Fig. [3j as a function of \(\ 2 for two typical values of the Majorana 
mass scale [Ib. = xn.mw- ^ n eacn case i we show the values for two choices of |mjv|/(U7?. = 10~ 4 , 10 -5 , given by the solid 
and dashed lines, respectively. The dotted horizontal line shows the experimental upper limit on the active neutrino 
mass, assuming a normal hierarchy with vanishing v e mass. From Fig. [3] we see that even in this simple scenario, it 
is possible to have a sterile neutrino in the keV-mass range to be a DM candidate (for \m^\j ' \ir ^ 10 -4 ) or in the 
eV-range to explain the LSND anomaly (for \itin\/hr & 10 -5 ). In both the cases, we must have |C| 2 -C 1, so that the 
lighter mass eigenstate is mostly the active neutrino (cf. Fig. [2]) with mass in the sub-eV range. 



IV. THE THREE FLAVOUR CASE 



We now turn our attention to the general case with three active neutrinos (z^ e ,i^,^ T ) and three light sterile 



neutrinos (y Sx , v S2 , , 
unitary matrix U$: 



In this case, we diagonalize the full 6x6 light neutrino mass matrix in ( 10 1 by a general 



M 



elf 



■Ut 



Mi 



2:s 



0, 



3 M 456 



6 ! 



(16) 



() 



where M123 
Moreover, 



diag(mi, m^) and M, 



150 



diag(w4, m5, mg), with mi,. ..,6 being the mass eigenvalues of A4 s- 



Ik 



( U el 

u Sll 

U S2 i 

\ Us 3 l 



U e 2 U e 3 U e 4 U e 5 

Up,2 Ufj,3 U^4 U^s 

U T 2 U T 3 U T i U T 5 

U Sl 2 U Sl 3 U Sl 4 U Sl 5 

U S2 2 U S2 3 U S2 4 U S2 5 

U S3 2 U S3 3 U S3 4 U S3 5 



U e6 \ 

u Sl6 

U S2 6 
U S3 6 J 



(17) 



is a 6 x 6 unitary matrix that relates the flavour states {u e , v^, v T , v Sl , v S2 , v S3 } to the mass eigenstates -{Vi, V2, ■ ■ ■ , v&}- 

As shown in Appendix A, Uq can be parametrized in terms of 15 Euler angles and 10 Dirac phases 3 . This, together 
with the six mass eigenvalues, introduces a total of 31 parameters on the right-hand-side of ( 16 ) some of which can 
be constrained experimentally, as follows, 
three mixing angles #12, #23, 



The only experimentally measured neutrino parameters so far are the 
?i3 between the mostly active neutrinos, and the two mass-squared differences Am^ 
and Am^ £Q. Also, there is some hint for a nonzero Dirac CP phase 1^13 from the global analysis of the 3- neutrino 
data 2, 3 . For our numerical analysis, we use the latest best- fit values of these parameters as given in |2], taking into 
account the reactor SBL data and assuming a normal hierarchy: 



'21 



33.3° , 623 - 
= 7.50 x 10" 



40.0° . 
' eV 2 , 



P13 
A ™31 



8.6° 



(pl3 = 1.677T, 



2.47 x lO' 6 eV z 



(18) 



On the other hand, for the mass and mixing parameters involving the sterile sector, there has been a number of 
recent analyses [I0HT4] involving one, two and three sterile neutrinos, denoted as (3+1), (3+2) and (3+3)-scenarios 
respectively. Although their global best-fit values for any particular scenario differ somewhat from each other, they 
all agree on the conclusion that there is a strong tension in the global data for the (3+l)-scenario [51123. This is 
mainly because the neutrino versus anti-neutrino and appearance versus disappearance data sets prefer two different 
values of Am| x for the (3+l)-scenario |13j . The neutrino versus anti- neutrino discrepancy is somewhat reduced in the 
(3+2)- and (3+3)-scenarios due to the additional CP-violating phases, thus leading to better global fits. However, the 
appearance versus disappearance tension still remains, and in particular, the fit to the latest MiniBooNE low-energy 
data is not improved by adding more sterile neutrinos. Thus due to the lack of a clear preference at the moment, we 
will consider benchmark cases for our model using both (3+1)- and (3+2)-fits, and will also comment on the (3+3)-fit. 

There also exist upper bounds on the absolute neutrino mass scale which could in principle be used to constrain the 
masses of the light sterile neutrinos and their mixing with the active ones. These constraints are briefly summarized 
below: 



1. The effective neutrino mass parameter in /3-decay experiments: 



lj2\Ue. 



(19) 



which determines the distortion of the electron energy spectrum due to nonzero neutrino mass and mixing. The 
current most stringent upper limit on this parameter is mp < 2.2 eV 30 31 . which is still compatible with the 
entire favoured region of the global fits [32] . The upcoming KATRIN experiment [33] with estimated sensitivity 
reach, as low as 0.2 eV, should be able to probe a substantial fraction of the allowed parameter space [51] . 

2. The effective neutrino mass parameter in neutrinoless double beta decay (0^/3/3) experiments: 



\U R . 



(20) 



3 For simplicity, wc set all Majorana phases to zero. 
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where a, are the Majorana phases associated with the elements U e i- The strongest bound so far on (m)pp < 
0.26 eV 35-38J and significant improvements are expected to take place in the near future [3S]. However, since 
the Majorana phases in ( pO) are still unknown, it is possible to have cancellations among them [ID], and hence, 
the constraint on may not apply for the sterile neutrino sector. For our benchmark model parameters 

discussed in Section [vj the 0vj3j3 contribution from the sterile sector is known to be negligible [41] . 

3. The sum of the neutrino masses, rn,, contributing to the total energy density in our Universe. This parameter 
can be constrained from the analysis of data on several cosmological observables [32], and the current upper 
limits arc in the eV range which are in tension with the sterile neutrino-favored region of the SBL data. The 
cosmologically-derived bounds, however, depend on the assumptions made on the cosmic history of the early 
Universe, as well as on the corresponding cosmological data set used, besides the uncertainties associated with 
the observables [TJ1 HFJ EH] • It might be possible to resolve some of these issues by analyzing the upcoming galaxy 
survey data (44) . It might also be plausible to consider a deviation from the standard cosmological picture, in 
order to accommodate the eV-scale sterile neutrinos as favoured by the SBL data [8]. A proper analysis of this 
possibility is beyond the scope of this paper. Hence, we do not consider the cosmological bounds for the rest of 
our discussion. 



In addition, the constraints on the active-sterile mixing from the non-unitarity of the PMNS mixing matrix [45] 
need to be taken into account as well. In the inverse seesaw model with /i^ = and fj,s ^ 0, the non-unitarity effects 
are proportional to ||£|| 2 [JSJ- For the MRISM with u fi ^ and [is = 0, the non-unitarity effects due to the mixing 
between the active and heavy sterile neutrinos in O) is of order ||Md/z^ || 2 which, as we will see in the following 
section, is very small (~ 1CP 9 ). This justifies our diagonalization procedure in (16) using a unitary matrix Uq. With 
the parametrization of Uq as given in Appen dix A , the non-unitarity due to the mixing between the active and light 
sterile neutrinos is determined by [cf. (A7)]. Using a combination of neutrino oscillation data and unitarity 



constraints in weak decays, the following allowed ranges for \V\ have been derived 



\V\ = 




0.65 
0.74 
0.75 



< 0.20 
0.57-0.82 
0.54-0.82 



(21) 



However, the constraints derived from weak decays are not directly applicable for the sterile neutrino masses below the 
electroweak scale as in our case. Also the constraints from Lepton Flavour Universality which may lead to measurable 
enhancements in various flavour physics observables [57] are not applicable to our case with all the light sterile masses 
much below the MeV scale. We find that the numerical values of the active neutrino mixing matrix elements in both 
of our benchmark scenarios are well within the la interval of the current global fit values from 3-neutrino oscillation 
data [2]: 



W\ 



0.795 - 0.846 
0.205 - 0.543 
0.215-0.548 



0.513-0.585 
0.416 - 0.730 
0.409 - 0.725 



0.126-0.178 
0.579 - 0.808 
0.567- 0.800 



(22) 



Note here that the constraints due to (22) are actually stronger than those obtained from (21) 



We conclude this section by deriving a relationship between the general unitary matrix that fully diagonalizes the 
effective neutrino mass matrix A^cff in ( 10 1 and the unitary matrix V introduced in (|6| which only block-diagonalizes 
A^eff as in (11). This relationship can be obtained by rewriting (A5| as follows: 



U 6 



1 3 Oc 

o 3 u. 



A R 
S B 



V 3 
o 3 1 3 



V 



i 3 o 3 

3 Ur 



V 3 
3 1 3 



(23) 



where V is identified as 



V = 



la 3 

3 U 



A R 
S B 



1 3 0; 

o 3 u, 



— I - L 3 



0, 



3 U 



V 



1 3 0; 

o 3 u, 



(24) 



and V has an exact representation as in (|6j with the replacement ( — !• £ = C^o- 



V. BENCHMARK SCENARIOS FOR LIGHT AND SUPERLIGHT STERILE NEUTRINOS 

In this section, we consider two physically interesting benchmark scenarios as shown in Fig. [4] and comment on 
a third one considered recently in the literature. We develop a method which enables us to fit some of the model 
parameters for each benchmark scenario, so as to be compatible with the SBL data. It also allows us to determine 
the values for the remaining theoretical parameters that have not yet been constrained by the existing neutrino data 
and could be probed in future experiments. 
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A. The (3+2+1) Case 



B. The (4+1+1) Case 



FIG. 4: The two benchmark scenarios discussed in Sections VA and VB The solid (dashed) horizontal lines represent the 



mostly active (sterile) neutrino mass eigenstates. The masses are not to scale. 



Employing the block form of IA§ in (A5), we obtain from (16) the following consistency conditions on the model 
parameters MjjjMjy and [Ir: 



R*M 123 & + U*M 4m U* = Mn^mJj, 
V*M 123 & +R*M i56 U^ = M D ^ l Ml, 
VMifflW + R*M i5e rf = M D (i^-(l 3 + x R f{x R j) M T D . 



(25) 
(26) 
(27) 
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FIG. 5: The magnitude of the mixing element (7 T 3 as a function of the mixing angle #34. The shaded region shows the 
experimentally allowed range. 



Using the best-fit values for the mass and mixing parameters discussed in Section [TV] we can derive constraints on the 
MRISM parameter space. For simplicity, we assume normal hierarchy 4 for the active neutrinos, i.e., mi < m-i JTI3, 
with mi = 0. 



A. The (3+2+1) Case 



In this scenario, we assume that all three active neutrinos are lighter than the three light sterile neutrinos, as shown 
schematically in Fig. [4]\. Among the three light sterile neutrinos, we assume one to be of order keV-scale to account 
for a warm/cold DM candidate and/or pulsar kicks [49]. Following the analysis of [50], we take the largest mass 
eigenvalue to be m 6 = 1 keV and the mixing parameters \U a ^\ = 7x 1CP 3 , for a = e, fi, r, which implies the mixing 



angles 9ie = 0.4°, for i = 1,2, 3 in ( A5 1. We also choose tp^ = for i = 1, 2, 3, 4 since these phases do not affect the 
SBL data significantly. The other two light sterile neutrinos are assumed to have eV-scale masses in order to explain 
the LSND+MiniBooNE results and reactor anomalies [10TU4| . Thus we have a simple extension of the (3+2)-models 
discussed in the literature which can now account for the DM as well. 

Since the third sterile state is much heavier than the other two and its mixing effects are very small, we can still 
use the global analysis of the (3+2)-scenario. From the latest global fit of [13] including the most recent MiniBooNE 
data 5 , we get the best-fit values of Am|i = 0.92 eV 2 , Amfi = 17 eV 2 , and 

|C7 e4 |=0.15, 1^1=0.13, |C7 e5 | =0.069, |^ 5 |=0.16, 5 4 = axg(UaU^U^U^) = 1.8tt , (28) 

which imply 

^15=4.0°, 2 5=9.2 O , 0i4 = 8.6°, 2 4 = 7.7°, </> 54 (<Pu ~ Pm) + (v>25 - Vis) = 1.8tt . (29) 

Note that the only experimentally relevant CP-phase in the (3+2)-case is the combination ^54, which enters the 
Vfj, — > D e transition probability. Thus, we cannot determine all the 5 Dirac phases in the (3+2)-sector. We assume 
for simplicity that the Dirac phases V24, Vis, V25, </?35 are all zero 6 so that ^54 ~ tpn. We also take the mixing 
angles $34 and $35 to be equal which can now be constrained from the allowed range of values for \U T 3\ since in the 



4 The case of inverted hierarchy can be worked out in a similar way. For a discussion of all possible scenarios with two and three sterile 
neutrinos, see |48| . 

5 The allowed region in the Am^ — Arrigj plane from the global fit of 1 131 is somewhat different from the other recent global fits |10f - 
1121 1 14| . This is mainly because the disappearance data sets used in 1131 prefer a medium Am^ (0.92 eV 2 ) and a high Am^ (18 eV 2 ) 
whereas the appearance data sets prefer a low Arre 2 ^ (0.31 eV 2 ) and a medium Am^ (1.0 eV 2 ). This incompatibility is expected to be 
addressed soon after more MiniBooNE neutrino data become available. 

6 We find that allowing a nonzero value for any of these phases does not change our fit significantly. 
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PMNS parametrization this element depends only on the mixing angles #23 and #13, whose values are now known 
experimentally [T . This is shown in Fig. [HJ where we plot the dependence of |C7 T 3| on # 34 — #35. The allowed range 
of \U T 3\ as given in (22) is shown as the shaded region. Thus, we see that in the first quadrant only 834 < 22.1° is 
allowed. 

While 6*34 and 645 are constrained to be small as shown in Fig. [5] the other three mixing angles, viz. #45, #46 and 
#56, parametrizing the mixing effects in the purely sterile sector are mostly insensitive to oscillation experiments. 
However, using the consistency conditions for our model parameters given by (25)— (271, we can predict their best- fit 
values along with the model parameters satisfying these conditions. 



Our approach to data fitting may be described as follows. For a given set of values for the mixing angles 
{#34; #35; #45; #46; #56} an d the scale /jr, we s olve for the mass matrix Mjv exactly using (25) for a symmetric 
structure of Mpr, 



i.e., for M/v = Mjy. Then using ( |26[ ), we obtain the input values for the elements of the mass matrix 
Md, which we use to check the validity of (27). More precisely we minimize the dimensionless function 



(y*M 123 v^ + R*M i56 R^ 


- 


M D ^(l 3 +x R f{x R )) Ml,] 




V*M 123 







(30) 



with respect to the set of variabl es { #34, #35, #45, #46, #56} while requiring that the elements U T i for i = 1,2,3 lie 
within the allowed range given in (22 ). For illustration, our results for x R = 1 case are shown in Fig. [6] and the best-fit 
values of the mixing angles that minimize the function J- given by ( 30 1 are 



#34 = 0.1°, #35 = 19.6°, #45 

which yields for the active neutrino mixing matrix 



26.4° 



'46 



12.9° 



54.8° , 



(31) 



0.8150 0.5354 0.1475 
|V| = ( 0.4793 0.5881 0.6182 
0.3032 0.5826 0.6783 



(32) 



and the model parameter values 



M D = 



M N = 



0.0605 - 





0210i 


0.0529 + 0.0105i 


0.0127- 


- 0.0066i 


0.0839 + 





0002i 


0.0743 


-0.0270 


- 0.0001^ 


0.0941 + 





0003i 


0.1321 - 0.0002i 


-0.1131 


+ O.OOOli 


0.7508 + 





0004i 


1.5916 + 0.000H 


0.9986- 


0.0003i \ 


1.5916 + 





OOOlz 


5.8432 - 0.0002i 


3.7990 + 


0.0002i 


0.9986 - 





0003i 


3.7990 + 0.0002i 


3.1780 - 


0.0002i / 



MeV, 



MeV 



(33) 
(34) 



with || CD = 0.44. Note that this leads to a very small non-unitary effect due to the mixing of the active neutrinos 
with the heavy sterile states: || A/d/i^ 1 1 ~ 10~ 9 . On the other hand, the non-unitarity due to the mixing with the 
light sterile states is well within the la uncertainty in the active neutrino mixing matrix, as can be seen by comparing 
J22| with (f32|. 



B. The (4+1+1) Case 



In this case, we assume that one of the sterile neutrinos is in the keV-range to account for the DM as in the 
previous case. Another sterile state is assumed to be in the eV-range to explain the SBL data. Thus, this can be 
regarded as a simple extension of the (3+l)-models discussed in the literature. The third sterile state is supcrlight 
(well below the eV range), having a weak mixing to the active sector, so that it does not contribute much to the 
SBL data. This scenario is shown schematically in Fig. [4}3. To the best of our knowledge, an explicit model with 
such a superlight sterile neutrino in the vicinity of the active neutrinos while being consistent with the SBL data has 
not been considered in the literature so far. Here, we would like to explore such a possibility and its experimental 
prospects. 
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FIG. 6: The function T denned in (30 1 versus the active-sterile neutrino mixing parameters for the (3+2+1) case. 



A superlight sterile state in the vicinity of the electron neutrino was considered in [ST] with Am^ ~ (0.2 — 2 



10~ J eV and the mixing angle sin 1Q\§ ~ 10 in order to explain the absence of an upturn at low energies in the 
solar neutrino data. They also required the mixing | C7^ s 3 3 1 2 ~ 0-02 — 0.2 in order to have a significant production 
of this superlight sterile state in the early Universe, which could generate the additional effective relativistic degree 
of freedom, AN e g ~ 1 as observed from the recent Wilkinson Microwave Anisotropy Probe (WMAP) data [TS] 7 , 



provided the mixing is large: |J7. 



s 3 3 I 



0.1 — 0.2. For illustration, let us consider a similar situation as in [51] . with 



Am 2 51 = 2 x 10~ 5 eV 2 . 



7 15 



= 025 = 2.0° 



while for the eV-scale sterile state, we use the best-fit value for the (3+l)-scenario [F3] : 



Ami 



0.92 eV 2 , \U e A = 0.15, \U, 



0.17. 



(35) 



(36) 



Also, we assume for simplicity that the Dirac phases </?24, if 15, </?35, ¥>i6> ¥>26 5 ¥>36i </>46 are ai l zero. The allowed 
range of values for \U T a\ will then constrain the values of the mixing angles #34, #35 and 63s to be small. For #34 = #35, 
this is shown in the left-panel of Fig. [7] where we have shown the allowed range of 636 — [0, 2tt/9]. On the other hand, 
if we choose #34 = #35 = #36, then the mixing angles are constrained to be < 22.1°, as in the previous scenario. This 
is shown in the right panel of Fig. [7] 



Using the consistency conditions for our model parameters given by (25|-(27) , we can derive the best-fit values 
for the set of six mixing angles {#34, #35, #36, #45, #46, #56}> as well as the corresponding model parameter values. 
To achieve this, we minimize the function defined in fl3 0|) w hile simultaneously requiring that the elements \U T i\ for 



i = 1,2,3 fall within the allowed range of | j given in 



In addition, we require that the element | t/s 3 3 1 should 



7 Note that in our case, the eV-scale sterile neutrino could also contribute to AiV e gf, depending on its thermal history at the Big Bang 
Nucleosynthesis epoch. 
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FIG. 7: The magnitude of the mixing element U T 3 as a function of the mixing angle #34 = #35 with #36 as a free parameter (left 
panel) and with 633 = #34 (right panel). The horizontal shaded region shows the experimentally allowed range. 



be within 0.1 — 0.2, as discussed above. For xr — 1, our results are shown in Fig. [8] and the best-fit values for the 
mixing angles which minimize the function T defined in ( 30 ) are 



#34 = 10.6° 
with |t/ S33 | 2 = 0.16 and 



= 10.. 



= 23i 



$45 — 82.3° 



\y\ 



0.8165 0.5364 0.1477 
0.4788 0.5948 0.6221 
0.2730 0.5498 0.6380 



74(i 



= 1.0° , 6» 56 = 50.8° 



(37) 



(38) 



The corresponding model parameters are 



M 



D 



M N = 



0.0075 + 0.0021z 0.2210 + 0.0154i 0.2239 - 0.0188i 
-0.0072 + 0.0002i 0.2129 0.2337 

0.0439 - 0.0002i 2.76678 - O.OOOli 2.3209 + 0.0002i 



0.0258 0.0910 0.1060 
0.0910 5.0303 3.8820 
0.1060 3.8820 3.4258 



MeV. 



MeV, 



(39) 



(40) 



with IICII = 1.45. Again in this case, we find that the non-unitarity effects are small, as can be seen by comparing 
p2l with ([38]). 



C. The (3+3) Case 



Finally, we would like to comment on the case where all the three sterile neutrinos are in the eV-range and heavier 
than the active ones. It was shown in |13j that this scenario has the highest compatibility between the neutrino and 
anti-neutrino data sets, but still has a poor compatibility between the appearance and disappearance data sets. They 
obtained the following best-fit values for the experimentally relevant 12 parameters in the (3+3)-case: 

Amii = 0.90 eV 2 , Am 2 51 = 17 eV 2 , Am^ = 22 eV 2 , 

|t/ M |=0.12, |[/ e4 |=0.Tl, 11^51=0.17, \U e5 \=0,U, 1^81=0.14, \U e6 \ = 0.11, 

054 = arg(U e5 U; 5 U: 4 U^) = 1.6tt , 64 = arg([/ e6 [/; 6 [/ e * 4 (7 Al 4) = 0.28^ , 65 = arg^f/^f/^) = 1.4tt . (41) 
For the parametrization of the unitary matrix discussed in Appendix A, this yields 

9 16 = 6.3°, 26 -8.1 o , 0i5-6.3°, 2 5-9.9°, 9 U = 6.6° , 9 24 = 7.1° , (42) 

054 — (<Pl4 - ^24) + (<£>25 - Vis) , 064 — {fli ~ ^24) + (^26 ~ Pie) , 065 - (fl5 ~ ^25) + (</?26 ~ Vw) • (43) 
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FIG. 8: The function T defined in (30 1 versus the active-sterile mixing parameters for the (4+1+1) case. 



Thus, we find that the three experimentally relevant Dirac CP phases for the (3+3) scenario are not independent 
of each other, and only two of them are independent, as already noted in |52j . Hence, we cannot obtain a solution 
for the phases in (43) satisfying the values given in (41), which were presumably obtained by assuming that they are 
independent parameters. Note that it is also very difficult to make the (3+3)-scenario with the best-fit Am 2 values 
given by ( 41 ) compatible with a standard cosmological model. 



VI. CONCLUSIONS 



We have discussed the possibility of having light and superlight sterile neutrinos in the recently proposed Minimal 
Radiative Inverse Seesaw Model. In the limit 3> ||Md||, ||Mjv||j after integrating out the heavy singlet states 

with masses of order ||//r||, we are left with three light sterile neutrinos along with the three usual active neutrinos. 
We have considered two benchmark scenarios for the mass hierarchy of these light sterile states. In the first scenario, 
denoted as the (3+2+l)-scenario, one of the light sterile states has mass in the keV-range and very small mixing 
with the active neutrinos to account for the Dark Matter in the Universe, whilst the other two sterile states are in 
the eV-range with non-zero mixing with the active states, as required to explain the LSND+MiniBooNE+reactor 
neutrino data. 

In this article, we have also discussed another potentially interesting scenario, denoted as the (4+l+l)-scenario, 
with one superlight sterile neutrino state almost degenerate in mass with the solar neutrinos, within the context of 
the Minimal Radiative Inverse Seesaw Model. This superlight sterile neutrino may mix weakly with both solar and 
atmospheric neutrinos and might give rise to observable effects in the current and future neutrino data. Such superlight 
sterile neutrinos may also provide an explanation for the extra radiation observed in the Universe. Moreover, this 
scenario is more compatible with the cosmological constraints on the sum of the neutrino masses. 

Finally, we note that both of these scenarios discussed here constitute complete extensions of the Standard Model 
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for neutrino masses and Dark Matter, while being consistent with all the existing data. In addition, the heavy singlet 
neutrinos in this model with a degenerate mass spectrum could be used to explain the observed matter-antimatter 
asymmetry in the Universe, through the mechanism of resonant leptogenesis. 
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Appendix A: Parametrization of a general unitary matrix 



Following [531 , we write a general n x n unitary matrix as 

Un = (w n _l )n W„_2,n • • ■ <^l,ri)(w„-2,n-l ■ • -Wl,n-l) ■ • ■ {^2,^1^)^1,2 , 

where, for j — i = 1, the rotation matrices Wjj are real: 

/ 1 \ 



UJ, 



1 



V 

and for j ' — i > 2, the rotation matrices are complex: 
/I 



1/ 



cos) 



• 
1 



• sin( 



sin%e 1 ^ ■ 



• cos 9i 



(for (j-i) = l), 



(for (j - i) > 2), 



1 J 



(Al) 



(A2) 



(A3) 



For our case of our interest here, i.e., for n = 6, ( Al ) becomes 

^6 = (^56^46^36^26^16) ( w 45 w 35 w 25 w 15) (^34^24^14) (W23W13) W12 



(A4) 
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Since W45 commutes with ujie for i = 1,2,3 [cf. (A2 A3|], we can rewrite (A4) as 

^6 = (W56W46W45)(W36W26W16W35W25W15W34W24W14) (W23W13W12 



1 3 3 

3 U 



A R 
S B 



V 3 
3 1 3 



V R 
R U 



(A5) 



Here the 3x3 unitary matrices Uq,Vq involve the mixing only within purely sterile and purely active neutrinos 
respectively, whereas the matrices R, S are responsible for the active-sterile mixing. For convenience, we have defined 



V = AVo, U = U Q B, R = UqSVq. 



(A6) 



where V is the new Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix for the active neutrinos. The analytical 
form of the matrices in (A5), in terms of the mixing angles and phases, has been presented in |54) . 



Finally, we note that the unitarity of U§ implies: 

VV* =AA^ = 1 3 -RR^, 
U ] U = B^B = 1 3 - R^R. 

Thus, the matrix R measures the non-unitarity of the PMNS matrix V. 
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